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Abstract 

The current knowledge about fluctuation induced long-ranged forces is sum- 
marized. Reference is made in particular to fluids near critical points, for 
which some new insight has been obtained recently. Where appropiate, results 
of analytic theory are compared with computer simulations and experiments. 

I. INTRODUCTION 

Forces between particles are governed by fields which themselves can be considered as 
composed of particles mediating the interaction by continuous exchange processes. Most 
prominent in the macroscopic world are electromagnetic fields and gravitational fields. For 
simplicity we specialize to electromagnetic forces here, but the line of argument scetched 
out in the following is valid in general. Macroscopic bodies exert electromagnetic forces on 
one another whenever they are charged, but if they are neutral all electromagnetic forces 
apparently vanish. H.B.G. Casimir JTJ] was the first to realize that this is not quite correct, 
because the electromagnetic field is fluctuating. These fluctuations may be due to quantum 
fluctuations at zero temperature in vacuum or due to thermal fluctuations in a cavity which 
is in contact with a heat bath. In any case the fluctuation spectrum, i.e., the energies 
which are associated with the Eigenmodes of the system and the form of the Eigenmodes 
themselves are manifestly influenced by the geometry of the system. The geometry is given 
by an arrangement of surfaces which impose boundary conditions on the fluctuating field 
and thus determine its mode spectrum. The free energy, which contains all information 
about the thermodynamic properties of the system, is essentially given by a sum over all 
modes and therefore the free energy will become geometry dependent. If, for example, two 
uncharged metallic bodies are placed at a certain distance in vacuum the free energy of the 
configuration depends on the shape of the bodies and the distance between them. Therefore, 
there will be an effective force acting between the bodies, which is given by the derivative of 
the free energy with respect to their distance. Note that this force is a direct consequence 
of the influence of the bodies on the electromagnetic fluctuation spectrum. Apart from the 
macroscopic length scales set by the geometry there are no other length scales in the system 
which limit the maximum wavelength of the fluctuations and therefore the force is governed 
by powers of the imposed length scales and scaling functions of their ratios, i.e., the resulting 
force is long - ranged. Due to the history of their discovery these forces are now known as 
Casimir forces and the influence of boundaries on the functional form of the free energy is 
known as the Casimir effect. It should also be noted that the presence of additional bodies 
in the above setup modifies the force between any two of them, i.e., it is not possible to 
express the Casimir effect as a sum of two - body contributions only. 
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A. A classic example 



There is a vast body of literature on the various aspects of the Casimir effect in electro- 
magnetism which are beyond the scope of this article. For a summary we refer to review 
articles dedicated to these subjects @ and for a recent survey in a more general context we 
refer to Chapter 3 of Ref. ||. However, for the understanding of the mechanism and the 
interpretaion of the results it is instructive to demonstrate some of the fundamental physical 
concepts by a simple example, which we call classic here for its historic meaning. 

Let us assume that we have two parallel perfectly conducting plates placed at a distance 
L in vacuum. We further assume that the plates have an infinite lateral extension so that we 
consider the thermodynamic limit with respect to the surface area A of the plates and only 
discuss the free energy T per unit area, i.e., we calculate F = lim^oo J 7 /A. Note, that L is 
the only macroscopic length scale in the problem. The mode spectrum of the electromagnetic 
field in this parallel plate geometry is particularly simple. In three dimensions the wave 
vector q = (p, k n ) consists of a lateral component p = (p x ,p y ) which is unconstrained by 
the geometry and a discrete perpendicular component k n — rnr/L for n = 1, 2, 3, . . . due to 
the condition that the electric field vector at each of the metallic surfaces must be aligned 
with the surface normal. A single mode is then characterized by p and n and its energy 
level spacing is given by 

£ p , n = frcyp 2 + (wr/L) 2 . (1.1) 

The energy content E p>n of a particular mode is given by its occupation number m p n = 
0, 1, 2, ... in the form E pn = e p>n (m Pin + . The free energy per unit area is then given by 
(see also Ref. ||] for a recent reconsideration of the Casimir effect at general temperature) 

F = J T^E £ p,n + 2k B T j -fJL £ In [1 - exp(-e p , n /(k B T))] , (1.2) 



n=l 



where k B and T denote the Boltzmann constant and the temperature, respectively. An 
additional factor of two coming from the summation over the polarizations has already been 
included in Eq. (|1.2j ). The integration over p is carried out to an ultraviolett cutoff A and 
the sum is truncated at some maximum mode number N. The ultraviolett cutoff parameter 
A is typically determined by the radius of the first Brillouin zone of the plate material. If a 
is the lattice constant of the material we identify A = tt/cl. The maximum mode number N 
can be written as N = L/b, where b is also a microscopic length scale (see below). 

For T = only the first term in Eq.( [1.2[ ) remains and this is the first example studied 
by Casimir [|IJ. Here, we will discuss the high temperature limit k B T ^> hcA, because this 
allows us to illustrate some aspects of the calculations involved within continuum models 
like, e.g., the Ginzburg - Landau model on a rather elementary level. For the more general 
case of layered dielectrica at finite temperature (dispersion forces) we refer to the classical 
literature [S,BJ and to Ref. M. To leading order in hcA/(k B T) we obtain 

F = 2k B T / /4 V In (1.3) 
J (2vr) 2 ^ k B T 1 J 
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Note that the integral and the summation in Eq.( |1.3| ) are only meaningful for finite cutoff 
parameters A and N. However, we only need the final result in the limits AL ^> 1 and 
N 3> 1. For simplicity we identify k B T ~ hcir/b to the order of magnitude which implies 
AL <C N, i.e., b <C a. The integral in Eq. ( |1.3|) is elementary and the resulting terms can be 
arranged as 
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The sum over the terms in the second line of Eq. ( p..4[ ) converges so that we can immediately 
perform the limit N — > oo here. We obtain 
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where terms which vanish in the limit A" — > oo have already been dropped. In order to 
evaluate Eq.( p..5| ) further for large A" and AL we employ Stirlings formula and the series 
expansion of the logarithm. With A = ir/a and N = L/b we obtain the final result 
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where ((3) ~ 1.202 is a special value of the Riemann zeta function and the dots indicate 
contributions which are exponentially small in AL. All terms which vanish in the limit 
N oo (6 — > 0) have consistently been dropped. 



The decomposition of the free energy per unit area given by Eq.(O) is a special case of 
the general decomposition 



F = LF b + F s , a + F Sjb + SF t 



ab 



;i-7) 



for a film with two surfaces of type a and b. The leading contribution to F is proportional 
to L and it corresponds to the bulk contribution of the free energy. In our example it is 
given by the radiation pressure F& between the plates. The second contribution to Eq. (|1.6|) 
is independent of L and it therefore corresponds to the sum of the surface free energies 
or surface tensions F s a + F s ^, where a = b in the above example. The third contribution 
varies as L~ 2 and it corresponds to the fluctuation induced long - ranged Casimir interaction 
between the plates, which is the most prominent contribution to the finite-size part SF a f, 
of the free energy in our example. Note that the Casimir contribution is independent of 
the microscopic cutoff parameters a and b. Its absolute strength at a given distance L and 
temperature T is governed by the numerical constant A = — £(3)/(87r) ~ —0.0478 which is 
usually called the Casimir amplitude. The Casimir amplitude is negative here, so that the 
Casimir force between the plates is attractive. The Casimir interaction can be obtained in 
very elegant ways known as zeta - function regularization, algebraic cutoff, or exponential 
cutoff schemes 0]. Their equivalence with respect to the cutoff - independent Casimir 
interaction has been explicitly shown for the example presented here ||. For further details 
see also Ref. PI. 
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B. Critical phenomena and correlated fluids 



The above example for the Casimir effect appears to be very specific at first sight, but 
the functional form of the free energy given by Eqs.( |1.3|) and (|1.6| ) is far more general 
than it seems. In fact, the underlying mechanism which leads to fluctuation induced long - 
ranged forces only requires a fluctuating field with geometric restrictions and a macroscopic 
length scale L imposed by the geometry as the only limiting factor for the wavelength of 
the fluctuations. Any system which is at a critical point also meets this requirement. The 
fluctuating field in this case is given by the order parameter and each of the individual 
contributions to the free energy as given by Eq. ( |1.7| ) is a sum of a regular part and singular 
part which contains the critical behavior of the system. Right at the critical point the 
correlation length £ is infinite so that the distance L between the system boundaries provides 
the only macroscopic length scale as required for the occurrence of long - ranged fluctuation 
induced forces. Above the critical point the correlation length is finite and therefore the ratio 
£/L governs the range of these forces. The existence of Casimir forces in critical systems 
was anticipated by Fisher and de Gennes []nj in the framework of the so-called distant 
wall corrections to critical profiles, where in many cases the Casimir amplitudes govern 
the leading distant-wall correction term to the profile in the vicinity of one of the system 
boundaries. For details and an extendend list of further references the reader is referred to 



chapter 4 of Ref . M and to Ref . [ I I 



It is important to realize that the Casimir amplitudes A and the associated scaling 
functions 9(L/£) that take the place of these amplitudes for finite L/£ (see Ref. and 
Sec. 2) are universal, i.e., they do not depend on microscopic details of the system under 
consideration. Note, however, that the precise form of the scaling functions 9 depends on the 
definition of the correlation length £. For systems with surfaces the concept of universality 
classes raises the question of whether there is surface critical behavior and to what extent 
it is governed by universal surface critical exponents. During the 1980 's this question was 
answered in favor of the general ideas of critical behavior and universality, i.e., microscopic 
surface properties are indeed unimportant. One only has to specify the type of boundary 
condition which the surface imposes on the order parameter. In this respect there are only 



three fundamentally different surface universality classes [12]. In particular, the surface may 



enhance the order parameter with the result that the system undergoes a second order phase 
transition in presence of an already ordered surface (extraordinary transition, E). The surface 
may also suppress the order parameter with the result that the system undergoes a second 
order phase transition in presence of a disordered surface (ordinary transition, O). Finally 
the surface and the bulk may order at the same temperature (special transition, SB), so the 
critical point of the system is in fact a multicritical point. If the spatial dimensionality d of 
the system is high enough there are two options for the occurrence of surface order above 
the bulk critical temperature T c ^. The surface may order spontaneously at a certain critical 
temperature T C)S > T Cjb or the surface may be ordered externally by the presence of a surface 
field. The bulk transition in presence of an externally ordered surface is called normal 
transition. However, it has been shown recently by rigorous arguments that the normal 
and the extraordinary transitions only differ by corrections to scaling, so both belong to 
the extraordinary surface universality class [[U|. Surface critical behavior has already been 
extensively reviewed 1^] (see also chapter 2 of Ref. for a short summary), so we refrain 
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from giving further details here. 

The distinction between the surface universality classes is vital for the Casimir forces, 
because the Casimir amplitudes A and the scaling functions 9 depend on them. The simplest 
boundary conditions apart from periodic ones are Dirichlet boundary conditions which sup- 
press the order parameter to zero at the surface. A system with these boundary conditions 
provides a representation of the ordinary surface universality class. The first systematic 
field-theoretic calculation of a Casimir amplitude was done by Symanzik [jTJJ] for this case. 
Starting from Eq. ( |1.3| ) the above example essentially reproduces all neccessary steps for such 
a calculation at the one - loop level (Gaussian theory). In general, concepts of the field - 
theoretic renormalization group are required which cannot be described here. The applica- 
tion of field theory to the critical behavior of finite systems is a field of ongoing research fl5| 
which has recently furnished unexpected results concerning the occurrence of nonuniversal 



critical finite - size behavior above the upper critical dimension [ 16[ . For reviews about the 



general concept of critical finite - size scaling the reader is referred to Ref. |T7[ . 

Above the critical temperature the range of the Casimir force is always limited by the 
correlation length, but below the critical temperature the situation may be different. In 
Ising - like systems the correlation length is also finite below T c ^ and therefore the Casimir 
forces have a finite range. If the system has a continuous symmetry, however, Goldstone 
modes cause correlation functions of the order parameter to remain long - ranged below 
T c fe. The most prominent examples are XY and Heisenberg ferromagnets which posess an 
O(N) symmetry with N = 2 and N = 3, respectively, in contrast to the Ising ferromagnet 
(N = 1). In other words, the correlation length of continuous ferromagnets remains infinite 
below T c b and therefore Goldstone modes also give rise to fluctuation induced long - ranged 
forces between system boundaries. Fluids with this property are sometimes called correlated 
fluids. The most important examples with respect to experimental realizations are liquid 4 He 
below the superfluid - normal transition |18| and nematic liquid crystals in the nematically 
ordered phase |nj, where fluctuations of the nematic director are responsible for the long 
- ranged nature of the Casimir force. Near the phase transition to the isotropic phase 
fluctuations of the degree of nematic order and the degree of biaxiality generate short - 
ranged corrections to the Casimir force |^0| . 

In summary, we have mentioned three options for the occurrence of fluctuation - induced 
long - ranged forces: the presence of long - ranged interactions (e.g., electromagnetism, see 
Sec.l.A), the presence of critical fluctuations, and the presence of Goldstone modes. In the 
following overview we will only consider the second option, i.e., systems in the vicinity of 
critical points. In particular, recent progress in the theoretical understanding of critical 
Casimir forces for all surface universality classes and especially for curved geometries will 
be presented. Special attention is also paid to the comparison of Casimir amplitudes and 
corresponding scaling functions with computer simulations and experiments, which are still 
in progress at this time. Due to the limited scope of this article other interesting develop- 
ments in related areas will not be described in any detail and an apology is made in advance 
to all authors whose work is not explicitly mentioned here. The remainder of this article is 
organized according to the three main approaches to critical Casimir forces, namely analytic 
theory, computer simulation, and experiments. 
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II. ANALYTIC THEORY 



The analytic theory of the Casimir effect in critical systems is based on the concept of 
finite - size scaling |3|,|17]]. Exact solutions of model systems in statistical mechanics give only 
limited access to the finite - size scaling functions, because they are mainly restricted to two 
dimensional systems. In d > 3 dimensions only the spherical model can be analyzed in a 



rigorous fashion j21[ which has recently been done with special regard to the film geometry 
for d = 3 dimensions f22| , |2"3"|| . Despite their limitations exact solutions provide valuable 
insight into the structure of the scaling functions and sometimes the results for d = 2 can 
be used to improve estimates obtained by approximative methods for d = 3 (see Sec. 3 of 
Ref. @ and below). 

The concept of finite - size scaling is a natural extension of the principle of scale invariance 
to critical systems with geometric constraints on macroscopic length scales. The principle 
of scale invariance itself may be viewed as a special case of the more general principle of 
conformal invariance (see Sec. 3 of Ref. || and Ref. |24j]). Conformal invariance implies the 
equivalence of systems with boundaries at T = T Cjfe if these systems can be mapped onto 
one another by a conformal transformation. The principle of conformal invariance holds in 
any dimension, but it is particularly powerful for d = 2 due to the exceptional large number 
of conformal mappings in this case (large conformal group, see Ref. [^|). Note that scale 
transformations are just very special conformal mappings. In the framework of conformal 
field theory the stress tensor plays a key role p3L|25|]. Here we only mention that the thermal 
average of the stress tensor yields the local Casimir force in a critical system and therefore 
the stress tensor provides a very important tool in the analytic theory of the Casimir effect. 
In fact, most of the Casimir amplitudes for d = 2 have been obtained from conformal field 
theory rather than exact solutions (see Sec. 3 of Ref. 



11) 



and Ref. 

Many of the experimentally relevant results have been obtained from a field - theoretic 
analysis of the well - known Ginzburg - Landau Hamiltonian 7i with geometric constraints 
which can be decomposed according to Ti. = Tib + 7~£ s + 7~£ e - The bulk contribution Tib is 
given by 



d d x 



I( V $) 2 + -$ 2 + ^ 
2 v ' 2 4! 



H ■ $ 



(2.1) 



for systems characterized by a iV-component order parameter $ = (</>i, . . . , <f>^) confined 
to a volume V, where N = 1,2,3 characterize the Ising, XY, and Heisenberg universality 
class, respectively. The parameters r and H correspond to the bare reduced temperature 
and external field. The physical (renormalized) reduced temperature and external field will 
be denoted by t and h in the following. The surface contribution 7i s can be written as 



d 



d-l 



X 



Hi ■ $ 



(2.2) 



where c and Hi correspond to the surface enhancement and the surface field, respectively 
T2|| . Note, that the surface 5* may consist of serveral disjoint parts. The last contribution 
7i e contains edge and curvature contributions to the Hamiltonian 7i which have first been 
considered in Ref. []rjj within the framework of the renormalization group. For experiments 
the ordinary transition (c = oo) and the extraordinary transition (e.g., c = — oo, see also 
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Ref. fL3"l ) are the most important cases. Neither TC S nor 7i e needs to be considered here in 
any detail, because their effect is completely contained in the boundary conditions for the 
order parameter $. We will therefore restrict the following discussion to the ordinary and 
the extraordinary surface universality class and to periodic boundary conditions. 



A. The spherical model 

The spherical model can be considered as the iV — > oo limit of O(N) symmetric classical 
spin models and it can also be expressed as the iV — ► oo limit of Eqs.( |2.1| ) and (|2.2| ). We 
only summarize the most recent results here, for a brief overview the reader is referred to 



Sec. 2. 2 of Ref. || and to Ref. pi| . In the presence of an external field h and for sufficiently 
small values of the reduced temperature t — (T — T C)b )/T c>b the singular contribution 5f ab to 
the finite-size part 5F ab of the free energy per unit area in a film geometry (see Eq.( [L.7p ) in 



d dimensions can be cast into the scaling form [22 



Sf ab (t, h, L) = k B T c ^L-^9 ab (tl}l\ hlPl v \ (2.3) 



near the critical point given by t = and h = 0, where ah indicates the combination of 
surface universality classes at the two surfaces. The critical exponents v and (3 characterize 
the temperature dependence of the correlation length £ ~ t~ v , t > and the order parameter 
(spontaneous magnetization) m ~ {—t)P, t < for h = 0, respectively. The exponent 5 
characterizes the functional dependence of the magnetization m ~ Ih] 1 ^ 5 on the external 
field h for t = 0. The form of the scaling arguments in Eq. ( |2.3j ) is imposed by the principle 
of scale invariance. They can be obtained by observing that L/£ is equivalent to the first 
scaling argument and L/£h, where £/, ~ h~ u ^ /36 " ) is the correlation length for finite field h at 
t = 0, is equivalent to the second scaling argument. For nearest neighbor interactions for 
d = 3 the critical exponents u, (3, and 5 of the spherical model are given by v — 1, /3 — 1/2, 
and 5 = 5. The special value A ab = 9 ab (0, 0) of the scaling function is the Casimir amplitude. 
In units of k B T cb the Casimir force K, ab = —-^Sf ab is characterized by the corresponding 
scaling form 

JC ab {t, h, L) = L~ d K ab {tl}'\ hL 135 ^) , (2.4) 
where the scaling function K ab is given by 

K ab (x, y) = (d- l)9ab(x, y) - -x-^-8 ab (x, y) - —y-^-6 ab (x, y). (2.5) 

v ox v ay 

The universal scaling functions 9 ab and K ab have been investigated recently for periodic 



boundary conditions ab = per [fffij23|| . For h = the scaling functions 6 per (x 1 0) and 
K per (x, 0) are both negative and increase monotonically with x, i.e., unlike the scaling func- 
tions in Ising - like systems they do not have a minimum in the vicinity of T = T C)b (x = 0) 
f22|| . For x — > +oo both scaling functions decay to zero exponentially, whereas for x — > — oo 
K per (x, 0) — > — C(3)/7T ~ —0.382. This behavior is due to the presence of Goldstone modes in 
the spherical model below T cb . For finite values of h (y ^ 0) the scaling functions per (x,y) 
and K per (x,y) again decay exponentially as L — > oo p2| . Moreover, both scaling functions 



also decay exponentially for y — >• oo at x = 0. The Casimir amplitude 
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Ap er $per(0, 0) 
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-0.15305... 



(2.6) 



can be obtained exactly p3[ and numerically it is very close the best availabe estimates for 
the Ising model for d = 3 (see Table | in Sec. 3). It has also been shown rigorously, that the 
scaling function 6 per (x, y) is a monotonically increasing function of each of its arguments as 
long as the temperature T is in the vicinity of T c ^ [23]. However, the hypothesis that this 



statement is true for any nearest neighbor O(N) symmetric spin model for N > 2 |23] could 
not be substantiated so far (see below). 



B. The Ginzburg Landau model 

1. Film geometry 

The film geometry has also been reinvestigated for the Ginzburg Landau model for the 



case of the extraordinary surface universality class [26], which is of particular interest for 
experiments with critical binary liquid mixtures. The scaling functions K^fa, 0) in zero 
external field have been determined within mean field theory for infinitely strong surface 
fields hi and h 2 which enclose an arbitrary angle a between a = (parallel surface fields) 
and a = it (antiparallel surface fields). For Ising-like systems only a = and a = ir can 
be realized and we refer to these cases as ab = ++ and ab = H — . The Casimir amplitude 
is negative for a = and positive for a = ir, it changes sign at a = n/3 [p6fl . Accordingly, 
the scaling function K ++ (x,0) is negative and the scaling function K + _(x, 0) is positive 
for all x within mean field theory, but it seems very likely that this is also true beyond 
mean field theory. The functional form of K ++ (x, 0) and K + _(x, 0) is illustrated in Fig.|l|, 
where the normalization of Ref. has been used. Note that both scaling functions take 
their extremal values at nonzero x, which makes them qualitatively very similar to the 
corresponding scaling functions for an Ising strip for d = 2, which can be solved exactly ]2"7| . 

The one-loop corrections to the mean field behavior are very hard to otain and at present 
they only exist for the Casimir amplitudes in the form of the ^-expansion, where e = 4 — d. 
The numerical quality of the ^-expansion when extrapolated to e = 1 is very poor so that 
exact results for d = 2 have been included in order to obtain an interpolation formula for 
the Casimir amplitudes for d = 3. The values for A ++ , A + _, and Aq+ obtained in this way 
agree reasonably well with Monte-Carlo estimates (see Sec. 3) fl2~6f . The e-expansion for 
A ++ and A + _ has recently received some independent confirmation from local-functional 
methods |28| which also provide reliable numerical estimates for d = 3 (see Table | in Sec. 
3). 

Apart from usual critical points, for which the upper critical dimension is d u = 4, tricriti- 
cal points in liquid mixtures with more than two components [25[] and in 3 He - 4 He mixtures 



(see Sec. 6 of Ref. 0) also provide an opportunity for experimental tests of the Casimir 
force. A theoretically appealing feature of a tricritical point is, that its upper critical di- 
mension is d u = 3 so that exact results for d = 3 can be obtained essentially from a mean 
field or a Gaussian theory. If, for example, Dirichlet boundary conditions are the correct 
ones for a 3 He - 4 He mixture in a film at the tricritical point, then the Casimir amplitude 
Aqo = — C(3)/(87r) ~ —0.0478 given in Eq.( |1.6|) is also the right one for this system. There 
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is, however, some debate concerning the correct boundary conditions for tricritical 3 He 
mixtures [29|. The result obtained for A ++ 



4 He 

3 contains a loga- 



at a tricritical point for d 
rithmic factor which is absent below the upper critical dimension d u and which introduces a 
dependence on a microcopic length scale into the Casimir amplitude |29|. This depencence 



is very weak and A ++ at tricriticality is expected to be about seven times larger than the 
corresponding amplitude at a usual critical point |29 . 



2. Curved geometries 

In view of experimental set-ups for, e.g., atomic force microscopy it is desirable to consider 
other geometries than films, because two plates cannot be kept parrallel accurately enough 
during force measurements. Curved geometries like a sphere in front of a planar wall or 
two spheres are much more convenient to control experimentally and are also much closer 
to reality in, e.g., colloidal suspensions |30| (see also Ref. [0). Some theoretical effort 
has therefore been made on the investigation of these curved geometries, where conformal 



invariance considerations have proved to be a very powerful tool at the critical point |32 
If T ab {r, Ri, Rq) denotes the free energy of a critical fluid in which two spheres with radii 
Ri and R2 at a center - to - center distance r > R\ + R2 are immersed, then the Casimir 
interaction bT ab takes the scaling form 



&F a6 (r, R u R 2 ) = F a b{r, Ri, R 2 ) ~ ^ ab {r = 00, R u R 2 ) = k B T C:b F ab (K) , (2.7) 

where ab denotes the combination of surface universality classes and k is the conformally 
invariant cross ratio 

k = (2R l R 2 )- 1 \r 2 - Rl- Rj\. (2.8) 

Note that the cases of two separate spheres in an unbounded critical medium and a single 
sphere inside a critical medium of spherical shape are conformally equivalent and are there- 
fore governed by the same universal scaling function F ab (n) f3~2fl . For large Ri and R 2 at 
fixed surface - to - surface distance D = r — R± — R 2 one obtains from the limit of parallel 
plates |H 

5F ab (r, R u R 2 ) = k B T C:b S d A ab [2 (D/R 1 + D / R 2 )]- {d - 1)/2 , (2.9) 

where Sd is the surface area of the unit sphere in d dimensions and A ab is the Casimir 
amplitude for parallel plates. In the opposite limit r ^> R\,R 2 the presence of the spheres 
can be taken into account by the small sphere expansion |32|] which yields 



Aj Af f R 1 R 2 \ Xi 
B, h \ r 2 



5F ab (r, R u R 2 ) = -k B T c i-^P^ ( ) * , (2.10) 



where ip = <j) is the order parameter if both a = b = E indicate the extraordinary surface 
universality class. In this case the scaling exponent is the scaling exponent of the order 
parameter — /3/u (c± 0.517 for the Ising model for d — 3). If only one of the surfaces is not 
characterized by the extraordinary surface universality class, the operator ip is given by the 
local energy density <fi 2 and x$ is the corresponding scaling exponent x$i = d — l/v (~ 1.41 
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for the Ising model for d — 3). The amplitudes Af and Af are the amplitudes of the critical 
profiles (V ; (^))to/2 = Af (2z)~ x ^ , s = a, b of the operator ip in a seminifinite system bounded 
by a planar surface of type s. The amplitude is the amplitude of the tptp correlation 
function in unbounded space. Although none of these amplitudes is universal individually, 
their combination in Eq. (|2.10|) is universal and its value for various surface types is exactly 
known for the Ising universality class for d = 2. In d = 4 — e dimensions estimates can be 



calculated from a renormalization group analysis of the Ginzburg Landau model | 3~2|| . Note, 
that Eqs.( [2.7|) , ( [2.9|) , and (|2.10|) only hold at the critical point. The Casimir interaction 
according to Eq. ( |2.1(J| ) is in fact very long ranged. For the extraordinary surface universality 
class it decays about as slowly as the Coulomb interaction. In all other cases the decay is 
faster, but it is still slower than the decay of, e.g., dipolar interactions. 

The full functional form of the scaling funcions F ++ (k), F +S b(k), and F +0 (k) 

has been calculated within mean field theory from the stress tensor in the concentric sphere 
geometry [Q. As for the case of parallel plates F ++ and F +S b are negative (attractive 

Casimir force), whereas F H and F + o are positive (repulsive Casimir force). The boundary 

conditions ab = 00, OSB, and SBSB have been treated within the Gaussian model, 
where Foo(k) = Fsb sb(k) < and Fo sb(k) > has been found. Although the analytic 
information from mean field or Gaussian theory is quite limited, the combination of these 
results with exact results for d = 2 yields fairly reliable estimates for ab = ++, H — , +0, 
and 00 within the Ising universality class in d = 3 fSlifl . Higher order calculations beyond 
mean field or the Gaussian approximation, respectively, for the concentric geometry are 
extremely demanding and results are not available. Finally, we note that the sphere - planar 
wall (SPW) geometry can also be obtained from the concentric geometry by a conformal 
mapping 

So far conformal invariance could be used to obtain the scaling functions of the Casimir 
interaction for various geometries with spherical surfaces. If the correlation length £ is 
finite, conformal invariance does no longer hold. Moreover, if all length scales £, r, R\, 
and i?2 are comparable, small sphere expansions cannot be made any longer and a new 
calculation is required for every geometry. In this case even mean field results can only be 
obtained numerically ||34|| . So far, this has only been done in detail for the SPW geometry 
(Ri = R, i?2 — ► oo, D = r — R\ — R2 = const.) with ++ boundary conditions and at arbitrary 
temperature near the critical point for Ising - like systems |[34| . We restrict the discussion 
to the case T > T c ^, where the correlation length £ + = governs the decay of the order 

parameter correlation function in real space. The Casimir force can be cast into the scaling 
form HI 

)C ++ (t,D,R) = ^±KX + (x + = D/i +1 y + = R/( + ) , (2.11) 

where a corresponding scaling function K^ + (x-,yJ) governs the scaling behavior of the 
Casimir force below T c ^. The scaling functions are obtained from the mean field evaluation 
of the stress tensor which requires the knowledge of the order parameter profile within mean 
field theory. The order parameter profile is obtained from a numerical solution of the Euler 
- Lagrange equation for Eq.( |2.1| ) in presence of parallel and infinite surface fields which 
dictate the boundary conditions. The functional form of Kl + (x + , y + ) is illustrated in Fig.||. 
As for the case of parallel plates, the Casimir force is attractive and takes its maximum 
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value above T c ^. The true position of the maximum is somewhat concealed in Fig.^| due to 
the normalization factor A 5 / 2 , which is required in order to absorb the divergence of Kf + 
for A = D/R — > 0. In this limit the Derjaguin approximation becomes valid, where the 
Casimir force is represented as an integral over parallel plate contributions. Each of these 
"parallel plates" in the d -dimensional SPW geometry is an infinitesimal annulus of width 
dp and radius p which is located on the surface of a paraboloid in order to approximate the 
sphere near the wall. The distance of one of these annuli from the wall is then given by 
L(p) = D + p 2 /(2R), where the integration is performed from p = to p = oo |34| . Note 
that this approximation is only valid for forces which decay sufficiently fast as L(p) — > oo. 
The amplitude of the Derjaguin approximation to the Casimir force at T = T C) & is indicated 
by the open circle in Fig.|2], where all solid lines meet. The dashed line correponds to the 
small sphere expansion to leading order, where y + = 1/5 has been used instead of the correct 
choice y + = 0, for which a factor A 2 is required for proper normalization [34] . 



The presence of small external fields can be used to drive, e.g., a critical binary liquid 
mixture slightly away from the critical concentration. Within the small sphere expansion 
the Casimir energy between two spheres (colloidal particles) turns out to be nonsymmetric 
with respect to deviations from the critical concentration so that the Casimir force is en- 
hanced when the concentration of the component preferentially adsorbed by the colloids is 
reduced [34]. This asymmetry is consistent with the asymmetry found experimentally in the 
flocculation phase diagram of colloidal suspensions . 



The concentric sphere geometry for ++ boundary conditions has also been considered 
at tricritical points [29fl, where the principle of conformal invariance can be used as well (see 
Ref . |I3 ] ) . The expressions for the Casimir forces in the different geometries are similar to the 
ones obtained for critical points [33j apart from the logarithmic dependence on a microscopic 
length scale. The scaling function of the Casimir force depends on the conformally invariant 
cross ratio given by Eq. (|2.S| ). In the range of distances D, where force measurements with 
the atomic force microscope appear to be feasable, both Casimir and van der Waals forces 
are essentially governed by the parallel plate limit of the curved geometries studied here 
|29fl . Corresponding investigations of the Casimir forces away from the tricritical point have 



apparently not been performed. 

Finally, we mention that a diluted polymer solution may also serve as a critical medium 
which mediates long - ranged forces between, e.g., colloidal particles p5] . Systematic in- 



vestigations, however, are still at an early stage and the description of these is beyond the 
scope of this article. 



III. COMPUTER SIMULATION 



Computer simulations of forces in liquid films have been performed in the past primarily 
with regard to the microscopic mechanisms of friction, adhesion, and lubrication, where 
both Monte 



Carlo 36 



and molecular dynamics methods [[37] have been used (see Ref. 
H for background material and more details). With regard to the Casimir force in critical 
or correlated fluids the situation is less satisfactory. The computational effort involved in 
such calculations is substantial and consequently only very few Monte - Carlo studies of 
the critical Casimir effect exist. Only rectangular geometries have been considered so far in 
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d = 3, because the currently available system sizes do not provide sufficient resolution to 
investigate curved geometries. 



A. Casimir amplitudes 

The first systematic attempt to measure the Casimir amplitudes of Ising and Potts 
models in a film geometry is based on a splitting procedure for lattice models at criticality 
3~8fl . The systems contain M d_1 x L lattice sites, where an aspect ratio of M/L = 6 turns out 



to be sufficient to approximate the film geometry. In the lateral directions always boundary 
conditions are always applied. A seam is introduced into the system Hamiltonian, which 
continuously weakens existing bonds and simulaneously establishes new bonds until the 
lattice is cut in two halfs of size M d ~ x xL/2. Histograms taken in the seam energy give access 



to the change of the free energy as a function of the seam strength |p8 |, which finally yields 
the total change of the free energy when the lattice is cut in two. For periodic boundary 
conditions this method yields the Casimir amplitude A per directly. For other boundary 
conditions the knowledge of A per is required as an input information [38]. The method 



works very well for d = 2 for critical Potts models with q = 2, 3, and 4 and has subsequently 
been applied to the Ising model for d = 3 with periodic boundary conditions [B3] and 



with surface fields [ 2"6fl . A summary of the currently available estimates for the Casimir 
amplitudes from various sources is displayed in Table | which includes older Migdal-Kadanoff 
estimates taken from Ref. |39[]. Apart from the well known numerical uncertainties regarding 
the extrapolation of the e-expansion and the Migdal-Kadanoff renormalization scheme the 
agreement between the estimates for each of the amplitudes is encouraging. Especially 

for A ++ and A H , where the e-expansion and the Migdal-Kadanoff scheme are particularly 

unreliable, the other estimates are fairly consistent. There are still some prospects to improve 
the Migdal-Kadanoff estimates also for these cases, but final results are not yet available 
40|[ . It would also be desirable to obtain estimates for the Casimir amplitudes from a 



field theoretic calculation for d = 3 directly, but attempts in this direction have not yet 
been made. The Monte - Carlo estimates for A ++ and A + _ are obtained by extrapolating 
the individual data to inifinite lattice size [p6 |. For A H this works rather well, but for 



A ++ substantial systematic uncertainties remain and additional data for larger systems are 



required to obtain a reliable extrapolation (see Figs. 4 and 5 in Ref. PEfl). At present local 



functional methods as set up in Ref. [28j seem to provide the most reliable estimates for A_ 



and A + _, because the dimensional dependence of these amplitudes appears to be captured 
rather well by the local free energy functional. Finally, we note that the Casimir amplitudes 
may also be accessible by exploring the order parameter distribution at the critical point 



41]. So far this method has only been used for fully finite cubic geometries, generalizations 



to other geometries have not yet been explored. 



B. Off-lattice models and wetting scenario 



A great drawback of the Monte - Carlo method introduced in Ref. |}8j is, that it cannot 
be generalized to temperatures T ^ T c ^. The method is based on the measurement of free 
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energy differences, which correspond to linear combinations of the scaling functions at dif- 
ferent scaling arguments for T 7^ T c ^. Data of extremely high accuracy would be required 
to disentangle the individual contributions to the measured free energy difference. An alter- 
native approach is to mimic the complete wetting scenario (see Ref. ||42|| ) near the critical 



end point of the demixing transition in a binary liquid mixture in a computer simulation 
§43|| . The order parameter in this case is the concentration of the mixture (N = 1, Ising 
universality class) rather than the density difference between liquid and gas, which is usu- 
ally taken as the order parameter near the liquid- vapor critical point. In fact, temperature 
and pressure are adjusted such that the mixture is in its vapor phase very close to liquid- 
vapor coexistence but far away from the liquid- vapor critical point. The interplay between 
the interparticle potential and the interaction between the particles and an external wall 
(substrate) may lead to the formation of a macroscopic liquid wetting layer of thickness L 



on the substrate at some temperature T w below the liquid- vapor critical point [}Eq| . The 
problem in the preparation of such a complete wetting layer for a binary mixture is to find 
a system, i.e., parameter values for a simulation, such that the critical end point of the 
demixing transitions is inside the complete wetting regime, where the macroscopic wetting 
layer remains stable (see Ref. for more background information on wetting transitions). 
The Casimir effect associated with the critical demixing transition can then be studied in a 
liquid layer of thickness L. The suggestion to probe the Casimir effect in complete wetting 
layers near critical end points was first made by Nightingale and Indekeu |44j and was later 
worked out in more detail, as the first estimates for the scaling functions of the Casimir force 
became available f|5] (see also Sec. 6 of Ref. ||). The main objective of such a simulation, 
however, is to obtain more insight into the Casimir effect in a more realistic off-lattice model 
with Lennard - Jones interactions, which would be the typical situation in an experiment 
4"3] |. Simulations have been performed for a symmetric binary liquid characterized by the 



Lennard - Jones interparticle potential 

Uij{r) = 4e 4i 



_ ,12 / a ..\e 



(3.1) 



for two paricle species i,j = 1,2, where the choices <jij = a for the hard core parameters 
and en = 622 = e, 612 = 0.7e for the well depth parameters have been made. Note that with 
these choices the system is invariant under the species exchange 1 2 and therefore the 
chemical potentials of both species have been set equal fi\ — fi2 — A* from the outset. The 
simulations are performed in a box of size P 2 x D, where periodic boundary conditions are 
applied in the x and y directions and two hard walls are specified in the z direction, one at 
z = and one at z = D. The wall at z = is characterized by the attractive potential 



V(z) 



2 /a\ 9 /a\ 3 
15 



(3.2) 



which acts equally on the particle species [|43f| . The interparticle interactions are truncated 
at R c = 2.5a, whereas no range cutoff is employed for V(z). The phase diagram of the 
system is spanned by the parameters fi/ksT, e/ksT, and e w /k B T. The system sizes used 
are P = 12.5a, 15a, and 17.5a, where D = 40a in all cases P3"| . The parameters are 
chosen such that i) a complete wetting layer forms on the substrate and ii) the endpoint 
of the critical demixing transitions is inside the complete wetting regime (see above). In 
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order to obtain a sufficiently thick wetting layer (L > 10a) the undersaturation of the 
vapor Sp/ ' p must be tuned to about 10~ 3 . The data aquisition is strongly hamperd by large 
fluctuations of the liquid vapor interface position (capillary waves) which also lead to a 
substantial slowing down of the algorithm. Furthermore, the data are strongly affected by 
lateral finite-size effects, because the capillary wave fluctuations also limit the lateral system 



sizes attainable within reasonable computational effort ||43|| . 

The equilibrium thickness L of the wetting layer (see also Sec. 4) minimizes the effective 
interface potential f42| , |4"3"|| as a function of the test layer thickness /, which is a variational 



parameter in the spirit of mean-field theory. The effective interface potential is given by 

uj(1) = l(pi - p v )5p + a w i + o iv + Su(l), (3.3) 

where pi and p v denote the liquid and vapor densities, respectively. The interfacial tensions 
a w i between the wall (z = 0) and the liquid and a[ v between the liquid and the vapor do 
not depend on /. The last term 8uj(T) contains the contributions of all interactions across 
the wetting layer and it depends on the boundary conditions. By construction there are 
no surface fields acting on the model liquid which break the 1^2 symmetry between the 
particle species, i.e., Hi = (see Eq. fl2.2|) for N = 1). However, the wall potental given by 
Eq.( |3.2|) acts like a negative surface enhancement c (see Eq.( [2.2|) ), which supports demixing 
near the surface, whereas the liquid-vapor interface acts as a free surface (c > 0) due to the 
internal 1 <-> 2 symmetry of the model liquid. The surface universality classes should thus 
be characterized by the combination {ah) = (0+). At T = T cep 5uo{l) is therefore used in 



the form j|3| 



M J) = £+*»T.,(^± + ^=), (3-4) 

where the Hamaker constant W ~ 2.5&BT ce? , governs the van der Waals contribution to 
the interaction potential and the Casimir amplitude Aq+ — 0.2 (see Table |) governs the 
Casimir interaction for a symmetric liquid mixture f43|. Note that positive Casimir ampli- 
tudes (repulsion) lead to a critical thickening of the wetting layer, whereas negative Casimir 
amplitudes (attraction) lead to a critical thinning of the wetting layer. The last term in 
Eq. fl3.4p is governed by the Casimir amplitude A per ~ —0.15. It provides an order-of- 
magnitude account of the aforementioned lateral finite-size effects which can be treated as a 
shift of the undersaturation 5 p. In the limit P 3> / (film geometry), the critical thickening 
of the wetting layer is given by [|3|] 

L c /L = (1 + k B T cep A 0+ /W) 1/3 , (3.5) 

where L c = L(T cep ) and L is the equilibrium layer thickness outside the critical regime. The 
measured film thickness L(T) versus temperature is shown in Fig.[3]. The critical thickening 
of the wetting layer L c / L — 1 for the largest system (lowest curve in Fig.^) is of the order 
of 3% which is in rough agreement with Eq. fl3.5p . The apparent reduction of the critical 
thickening with incrasing lateral system size can be explained semiquantitative^ by the 
lateral finite-size correction included in Eq. ( |3.4| ). Further studies of off-lattice models like 
this are certainly desirable, however, algorithmic improvements for the treatment of capillary 
waves will be indispensable for future progress. 
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C. Lattice stress tensor 



The computer simulation of the complete wetting scenario is quite successful, but by 
its design it only gives an indirect account of the scaling functions K of the Casimir force. 
As already described in Sec. 2 the most direct access to the Casimir force is given by the 
thermal average of the stress tensor and it would therefore be most convenient to have a 
lattice expression for the stress tensor available for spin models. Such expressions can indeed 
be obtained and successfully used in Monte - Carlo simulations for lattice models for d = 2 
f46| . The basic idea behind the construction of the stress tensor is the same is in continuum 



theory: one calculates the response of the free energy to a nonconformal mapping of the 
system, e.g., an anisotropic rescaling of the coupling constants. For example, one may choose 
Jx = Jx(ty, J y = JyW with J x (0) = J y (0) = J and J x (\) = J y (—X) on a square lattice, 
such that the critical temperature does not depend on A. For the Ising model for d = 2 this 
procedure yields, e.g., the xx component of the lattice stress tensor at lattice site in 
the form fl6| 

txx(i,j) = (SijSt+ij - SijSij+x) , (3.6) 

where J' x (0) is the derivative of J X (X) at the isotropic point A = 0. The thermal aver- 
ages (t xx ) of Eq. (|3lf) and (T xx ) of the stress tensor in conformal field theory are related by 
(t xx ) = ci(T xx ) up to corrections to scaling, where a is exactly known for the 2d Ising model. 
For periodic boundary conditions in a strip geometry Eq.( |3.6|) can be used to measure the 
Casimir amplitude A per (i.e., the conformal anomaly number c) and two or more scaling 
dimensions for the Ising model and other models if Eq. ( |3.6| ) is generalized accordingly [f46 



For lattice models for d = 2 conformal field theory provides sufficient background infor- 
mation so that the desired quantities can be extracted from elaborate fit procedures [06 



Although Eq. ( |3.6|) can be readily generalized to d — 3, additional information from confor- 
mal field theory, which is vital for the data interpretation for d — 2, is no longer available. 
Furthermore, (t xx ) still contains surface contributions for nonperiodic boundary conditions, 
because the surface tensions will depend on the anisotropy parameter A even if the critical 
point does not. However, for periodic boundary conditions (txx) is at least proportional to 
the Casimir force and some preliminary studies for the XY model in an M 2 x L geometry 
for d = 3 dimensions look promising [j|7| , although high statistics is needed already for small 
systems. 



IV. EXPERIMENTS 

Experimental verifications of the Casimir effect in critical liquids are exceedingly dif- 
ficult, because data of high accuracy are required and both samples and apparatus must 
be prepared with great care. At present, two lines of approach are considered, namely the 
wetting scenario scetched already in Sec. 3 and direct force measurements by atomic force 
microscopes (AFM). 
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A. Wetting experiments 



For a wetting experiment in the vicinity of a critical point a fluid is required which 
posesses a critical end point on the liquid vapor coexistence line. One option for this setup 



is provided by He near its lower A point [48|. For this system the interaction part of the 



effective interface potential (see Eqs.(|3.3j) and fl3.4Q ) must be modified according to the 
universality class of the A-transition in 4 He (XY, N = 2). The boundary conditions at the 
two interfaces of the wetting layer seem to be very well approximated by Dirichlet boundary 
conditions (O surface universality class). This leads to |]45| J 



where L x ~ 193 A denotes the crossover length to retardation ||49|| , and 9qo is the scaling 



function of the Casimir potential for the ordinary surface universality class (see Eq.( [2.3| )). 
Note that W and L x depend on the dielectric properties of the adsorbate and the substrate. 
From Eqs.( |3.3| ) and ( |4.1[ ) one expects a critical thinning of the wetting layer thickness, 
because 9oo < 0. Note that the second scaling argument of O o captures off-coexsistence 
effects due to the undersaturation S/jl of the 4 He vapor. At the A point (t — 0, 8fi — 0) one 
has #oo (0,0) = ^oo — —0.024 which results in a critical thinning of ~ 0.3% for standard 
substrates like, e.g., copper [j45|j3"l|. In the experimental setup a stack of five copper capacitors 
is placed inside a cell which contains liquid 4 He at the bottom. The surfaces of the capacitor 
provide the substrate potential (see Eq. (|3.2] )) and their elevation h in the gravitational 
field controls the undersaturation 5/x ~ p v gh of the 4 He vapor. The layer thickness is 
obtained from high precision measurements of the capacitance of each of the capacitors as a 
function of temperature. The wetting behavior of 4 He is extremely sensitive to the surface 
morphology of the copper plates. In particular, microscopic scratches and dust particles 
lead to localized condensation of 4 He on the surface which results in an overestimation of 
the thickness. Moreover, surface roughness leads to an enhanced surface area which also 
increases the amount of liquid 4 He on the substrate. Even with the most advanced polishing 
techniques these effects cannot be avoided completely and therefore also the experimental 
verification of the DLP theory of dispersion forces || remains a challenge. Nevertheless, the 
experimental data of the film thinning show a pronounced minimum well below T\ which 
is given by the specific value x m = —9.2 ± 0.2 of the first scaling argument x = tL x l v 
in Eq.( |4.1| ). This value coincides with the minimum of the scaling function Koo{ x iU) of 
the Casimir force with respect to x. The experimental estimate of d(x) = Koo{ x iU — 0) 
extracted from the data is displayed in Fig.f|, which does not show the expected data collapse 
for the scaling function Koo- On the contrary, a systematic trend in the data as function 
of the elevations h of the capacitors is visible as shown in the inset of Fig.(|. One possible 
explanation may be given by off-coexistence effects, which would require the second scaling 
variable y = 5fi L^ 8 I V in Eq. (|4.1|) for data collapse. If a linear dependence of Kqo on V 
is assumed, the deviations from data collapse are indeed drastically reduced |[48| . Another 
option is provided by the introduction of a roughness correction factor as suggested in Ref. 
| |18| |, which leads to a similar improvement [|48|| . However, it is evident from Fig.^, that there 
are no appreciable deviations from data collapse for x > 0, where a quantitative prediction 
for Kqo{ x i V — 0) exists [f|^,||. The comparison is displayed in Fig.|^ which shows reasonable 
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agreement between the data and the prediction. Finally, we note that a finite thinning effect 
remains visible for temperatures further below T\, as one would expect from the presence of 
Goldstone modes fl~8| , |2"3f| . The overall shape of Kqo is manifestly nonmonotonic in contrast 



to the recently stated monotonicity hypothesis for 0(N > 1) symmetric spin models [23 



As second option for a wetting experiment in the vicinity of a critical end point is 
provided by binary liquid mixtures near the critical end point T cep (see Sec. 3) of the line of 
second order demixing transitions |50J (see Refs. [i2~,4~5 for complete phase diagrams). The 
physical situation is very much like 4 He near the lower A-point, except that both the bulk 
and surface universality classes are different here. The second order demixing transition 
is characterized by a scalar order parameter (concentration), so the system is in the Ising 
(N = 1) universality class. The substrate material as well as the liquid - vapor interface, 
which provide the boundaries of the system, usually show some preferential affinity for one 
of the two components of the mixture, i.e., the concentration (order parameter) departs 
from its critical bulk value in the vicinity of the surfaces. This phenomenon is known 
as critical adsorption (see, e.g., Ref. [12|) and it is captured by the extraordinary surface 



universality class. In the experiment []50j] a molecularly smooth (100) Si wafer (n-type, 
phosphorous doping) covered with a Si02 layer of ~ 2. Oram thickness is used, which is 
suspended vertically inside a pyrex cell. The elevation h of the substrate above the bulk 
liquid, at which the wetting layer thickness is measured, controls the undersaturation of the 
vapor. The reduction of temperature gradients for wetting agents other than superfluid 4 He 
is a quite demanding task and it substantially complicates the preparation of the cell and 
the sample ||50|| . In this experiment two organic mixtures have been used, namely methanol 
+ hexane (MH) and 2-methoxy-ethanol + methylcyclohexane (MM). In MH the methanol 
component is adsorbed at the Si wafer, whereas hexane is adsorbed at the liquid - vapor 
interface of the wetting layer. In MM the situation is similar: the 2-methoxy-ethanol is 
adsorbed at the Si wafer, whereas the methylcyclohexane is adsorbed at the liquid - vapor 
interface. The wetting layers of both mixtures are therefore characterized by the scaling 
functions of the Casimir potential. The interaction contribution 8u(l) to the effective 
interface potential u(l) in this case is assumed to be of the form ]50| 



MO = y~ Ae ~ l/d + ^f^+- (* i 1/v , o) , (4.2) 

where retardation and off - coexistence effects are neglected. The exponential contribution 
to Eq. (|4.2|) is due to the presence of the hard wall, which structures the adsorbed fluid over 
a molecular distance d. The critical temperature T cep is about 300f^ [50]. At the moment 



only mean field results ^6] and exact results for d = 2 J27] exist for the scaling function 
of the Casimir force. In order to obtain reasonable estimates also for d = 3 at least the one 
- loop corrections are required which only exist for A + _ = 2i^ + _(0,0) at the moment (see 
Ref. |2(J and Table |). From Eq.( [3.5|) and typical values for T cep and the Hamaker constant 
W one expects a critical thickening L c /Lq > 2 of the wetting layer, when the estimate 
A + _ ~ 2.4 (see Table D is used. As T cep is approached from above a critical thickening 
of the wetting layer consistent with this expectation has been found in the experiment and 
the data for K + _ (x, 0) are indeed consistent with scaling [|50j . As function of the scaling 
variable y = L/^ + the scaling function •§ + ^{y) = K + _(x = (£o~?/) 1/V , 0) is shown in Fig|| 
A comparison between i9 + _(y) at and away from the critical composition is shown in the 
inset for MH. The shape of shown in Fig.|6| resembles that of the mean field estimate 
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K + - {x, 0) for x > shown in Fig.|T]. However, the Hamaker constant W and therefore also 
the Casimir amplitude A + _, which have been extracted from the data, are much smaller 
than anticipated. The reason for this discrepancy has not yet been fully understood. One 
possible explanation could be that the parameters of the system are not in the complete 
wetting regime, i.e., only partial rather than complete wetting ]22] is achieved. Further 
studies are currently under way. 



B. AFM measurements 

As already mentioned in Sec. 3 the theoretical investigation of critical fluids in curved 
geometries is, inter alia, motivated by the prospects of measuring the Casimir force directly 



with an AFM. At the moment only exploratory results are available which have been 



obtained for the SPW geometry in liquid crystals (see Refs. | 19|j20| ). In this study a tem- 
perature controlled AFM has been used to measure the force between a sphere mounted on 
the cantilever tip of the AFM and a planar wall immersed in an 8CB liquid crystal near the 
isotropic - nematic phase transition. Above the transition in the isotropic phase an attrac- 
tive force of the order of 10 -10 iV at a distance D = lnm±0.1nm between the surface of the 
sphere (R = 5/xm) and the wall is detected only when the two surfaces are moving apart. 
This phenomenon is similar to the capillary force in AFM microscopy and it is interpreted 
as the adsorption of a nematically ordered layer of the liquid crystal on the surface of the 
sensing probe [|5T|. Slightly above the transition to the nematically ordered phase an addi- 
tional attractive force of the order of 10~ n iV is detected when the surfaces are approaching 
one another. This additional force is conjectured to be the Casimir force mediated by the 



onset of Goldstone modes of the nematic director field in the ordered phase |19],[5T], where 



the boundary conditions are supplied by the type of anchoring of the nematic director on the 



surfaces [19]]. Further quantitative studies of Casimir forces in critical and correlated liquids 
with this apparatus are certainly desirable. Finally, we note that the radiation pressure on a 
dielectric sphere in the evanescent field of totally reflected light has recently been measured 
using such AFM techniques [j52|| . 



V. PROSPECTS FOR FURTHER INVESTIGATIONS 

The theoretical knowledge about Casimir forces in critical and correlated fluids which 
has been accumulated during the last 10 years has become so detailed, that the stage is set 
for experimental tests of various kinds. Wetting experiments near critical end points have 
already proved to be a powerful tool to accomplish this goal for quite a variety of fluids. 
Further studies in this direction are certainly highly desirable and the prospects for them 
are very good despite substantial experimental challenges one has to face. From the existing 
theoretical work on curved geometries it has also become clear, that the Casimir forces in 
critical and correlated fluids are within reach of current AFM designs. Preparation of the 
samples and temperature stabilization of the sample and the instrument again pose major 
challenges for AFM force measurements, however, the prospects of probing the Casimir effect 
quantitatively are also very good. 
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Conversely, the experimental approaches to the Casimir effect also pose new theoretical 
challenges. The problem of substrate roughness in wetting experiments has already been 
mentioned above and for the case of quenched roughness theoretical results already exist WE 



However, one of the boundaries of a wetting layer is a free liquid - vapor interface, which may 
undergo large scale fluctuations due to capillary waves. What kind of corrections capillary 
waves as additional degrees of freedom impose on the critical Casimir potential is an open 
question. Quantitative estimates of these corrections are not only important for experiments, 
they would also aid the data interpretation of computer simulations for critical wetting 
layers. To what extent off-coexistence effects influence experimental and numerical wetting 
layer data is also a largely open problem. In this respect the recently explored numerical 
access to the stress tensor of lattice models may prove particularly useful |46|j47| . Finally, 
it should be mentioned that improvements of existing theoretical or numerical estimates 
for the scaling function K a ^ of the Casimir force in particular for the extraordinary surface 
universality class in various geometries are still needed in order to extract the Casimir effect 
from experimental data as reliably as possible. 

Although the history of the Casimir effect goes back more than half a century it has 
remained an active field of research. This article can therefore only provide a snapshot of 
current knowledge in this area rather than a completed picture. If this presentation could 
finally help to trigger or direct new research work in this field, then its main purpose would 
be fulfilled. 
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FIGURES 



FIG. 1. Scaling functions K ++ (x, 0) (solid line) and (x, 0) (dashed line) taken from Fig.l of 

Ref. (2(|. The x range influenced by the bulk critical point x = is very broad and the asymptotic 

decay for x — ► ±oo is dominated by an exponential. Note that K ++ (x, 0) and (x, 0) take their 

extreme values at x ~ 10 and x ~ —25, respectively. 

FIG. 2. Scaling function K+ + {x + ,y + ) as function of x + for various values of y + (solid lines). 
The prefactor A 5 / 2 absorbs the divergence of the scaling function in the limit A — > 0, where 
the Derjaguin approximation becomes valid (see main text). The fixed point value u* of the 
renormalized coupling constant is required as an additional normalization due to the mean field 
character of the calculation. The dashed line corresponds to the small sphere expansion, which is 
shown here for y + = 1/5. The exponential decay of the scaling function sets in at x+ ~ 4. 

FIG. 3. Thickness of the wetting layer as a function of temperature in reduced Lennard - Jones 
units along a path parallel to the liquid vapor coexistence line. Data are shown for each of the three 
system sizes studied. The results were obtained from multihistogram extrapolation of simulation 
data accumulated at three points on this path, corresponding to temperatures T = 0.946, 0.958, 
0.97 (see Ref. @). 

FIG. 4. Scaling function = Koo( x ,Q) as a function of x. The magnitude of the minimum 
increases systematically with the height h of the capacitor. The measured layer thickness L is 
roughly between 300A and 600A depending on the capacitor index 1-5. The inset shows the 
value of i?(a?) at the minimum vs. height. The uncertainty in the vertical scale is 2 - 10% (taken 
from Ref. g§). 

FIG. 5. Blow-up of the region x > in Fig|4]. Every other data point is shown. The solid line 
shows the prediction from Fig. 9 in Ref. |^5[ (taken from Ref. p8f). 

FIG. 6. Universal scaling function y = L/^ + for the critical Casimir force. The 

symbols represent data at fixed elevations h = 1.5mm (diamonds), h = 3.3mm (squares) for MM 
and h = 3. 4mm (triangles), h = 6.3mm (inverted triangles) for MH. In the inset the experimental 
•&+—(y) for the system MH with critical composition (squares), 5% excess hexane (circles) and 
10% excess hexane (triangles) at two different heights (3.5 mm (open symbols) and 6.0 mm (solid 
symbols)) on a silicon wafer is shown. 
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TABLES 



TABLE I. Casimir amplitudes for the Ising universality class for d = 3. The values labelled 
e = 1 are obtained by extrapolating the e-expansion for N = 1 to s = 1 p§. The values labelled 
d = 3 are obtained from Pade type approximants for d = 3 (e = 1) [^6|. The Monte-Carlo estimates 
obtained from the algorithm presented in Ref. [38| are labelled by 'MC. Statistical errors (one 
standard deviation) are indicated by the figures in parenthesis. The last two lines show estimates 
taken from Refs. p| and 28]. 
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